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Abstract: We assume the DE state equations w(a) = wo + w a {a P — a), and study the 
dependence of the constraints on wq and w a coefficients on the pivoting redshift 1 + z p = 
l/a p . Coefficients are fitted to data including WMAP7, SNIa (Union 2.1), BAO's (including 
WiggleZ and SDSS results) and Hq constraints. The fitting algorithm is CosmoMC. We 
find specific differences between the cases when i/-mass is allowed or disregarded. More in 
detail: (i) The z p value yielding uncorrelated constraints on wq and w a is different in the 
two cases, holding 0.25 and 0.35, respectively, (ii) If we consider the intervals allowed to 
wo, we find that they shift when z p increases, in opposite directions for vanishing or allowed 
i/-mass. This leads to no overlap between \—a intervals already at z p >~ 0.4. (iii) The 
WQ-w a constraints found by using any pivot z p can be translated into constraints holding 
at a specific z p value (0 or the z p where errors are uncorrelated). When we do so, most 
error ellipsoids exhibit a nice overlap, (iv) the only exception being the constraints found 
for vanishing v— mass by using z p = 0, which turn out to be significantly more restrictive 
than others (up to ~ 30%). 
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1. Introduction 

Owing to the conceptual problems of ACDM, a number of options for Dark Energy (DE) 
nature have been considered. In particular, DE could be a scalar field, necessarily self- 
interacting and possibly interacting with Dark Matter [1], or just a phenomenological 
consequence of large scale GR violations [2]. But neither these options, nor still more 
exhotic hypotheses [3], led to appreciable improvements of the fit between theory and 
data [4]. 

The problem has then been tackled from the phenomenological side, by testing whether 
any linear w(a), different from w(a) = —1, improves data fits. A possible option amounts 
then to express the linear laws through the equations 

w(a) =w + w a (1 -a) , (1.1) 

aiming then at testing how various sets of data yield constraints on wq and w a . Here o 
is the scale factor, normalized to unity at the present time. The same linear laws can be 
expressed also through the equations 

w(a) = w 0<ap + w a ,a p {a p - a) (1.2) 

which differ from (1.1) for selecting a non-vanishing pivoting redshift 

Z p = 1/dp - 1 , (1.3) 

while we put an extra index to the linear coefficients WQ t a p , w ajap to put in evidence that, 
when changing z p , their values change. The straight lines defined by eq. (1.1) and eq. (1.2) 
are however the same: any equation (1.2) turns into an equation (1.1) if we set 

w ,a p = w - w a (a p - 1) (1.4) 
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w W Q 
Figure 1: 2D marginalized likelihood, when taking z p = 0, on the plane w — w a at 65% and 95% of 
confidence for different data sets, when assuming massless (l.h.s.) or massive (r.h.s.) neutrinos. The 
coordinate scale being the same on both sides allows us to appreciate how the WQ-w a uncertainty 
increases when a j/-mass degree of freedom is considered. 

and w aAp = w a . Notice that this last identity does not imply that limits on w aAp are 
indipendent from the pivoting redshift. In the sequel, whenever this causes no confusion, we 
shall however follow the common use and call Wo, w a the two parameters in any expression 
(1.2). 

Linear laws can be fitted to data by using different a p values. Here we aim at testing, 
first of all, how compatible are results obtained when varying the pivoting redshift. 

We shall do so in two cases: either neglecting or allowing the option that M v = 
Ylv m v 7^ (the sum is extended to the mass eigenvalues for 3 standard neutrino flavors). 
Let us also remind that the neutrino density parameter 

n u h 2 = 1.08 x 10- 2 (A^/eV)(T 07 /2.73 K) 3 , (1.5) 

so that, when the dark matter reduced density parameter u c = U c h 2 is assigned, the 
neutrino fraction /„ = fi„/f2 c immediately follows. 

The value of a p can be selected so to have uncorrelated phenomenological density 
constraints on wq an( A w a . Here we also wish to put in evidence that: (i) the pivoting 
redshift yielding uncorrelated constraints is different, if f u = or can be ^ 0; (ii) also the 
dependence on a p of the wq interval compatible with data depends on the above option. 

2. Results for z p = 0. 

Let us then report, first of all, the results of Monte Carlo fits of DE state equations vs. data, 
performed by using the algorithm CosmoMC 1 [5] (May 2010 version); the CosmoMC code 
was integrated with the first version of the PPF module 2 [6] for CAMB 3 [11]. Fits were 
performed in respect to the following parameters: wq, w a (in eq. 1.2) and uj^ = f^/i 2 , 

1 http : / / www. cosmologist . info/cosmomc 

2 camb.info/ppf 

3 http://www.camb.info/ 
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Neutrinos 


Massive 


Neutrinos 




Data Set 


Wq rt la rfc 2(7 


W a ± 1(7 ± 2(T 


Wo i 1(7 ± 2(7 


W a ± 1(7 ± 2(7 


f v ± la ± 2(7 


CMB+SN 

....+BAO 

....+HST 


-0 qq+O.ll+0.41 
u - yy -0. 13-0.36 
i -, n +0.08+0.34 
i - lu -0.10-0.30 
i i q+0.08+0.34 
i - lo -0.10-0.30 


n nzL +0.55+1.25 
u - u *-0.38-1.69 
n oi +0.39+0.93 
u - o± -0.23-1.29 
n qi +0.40+1.00 
u - o± -0.25-1.32 


i nfi +0. 13+0.46 
1 - uu -0. 14-0.45 
i n q+0. 10+0.39 
1 - uo -0. 12-0.35 
i n7 +0. 10+0.38 
1 - UI -0. 11-0. 34 


n oc+0.68+1.56 
u - zd -0.44-2.14 
n q 7 +0.60+1.41 
u -°' -0.36-2.01 
n qi +0.56+1.37 
U - Oi -0.34-1.89 


n n/i^+°- 013 + - 052 

U.UU_ Q45-0 .045 

n n/in+°- 011 + - 043 

U.U1U_ Q40-0 .040 

n n-3fi+°- 009 + - 041 

U.UOU_ 036-0 .036 



Table 1: Mean, la limits and 2a limits on Wq, w a and /„ (when relevant), with and without 
massive neutrinos and different data sets. Values are obtained by fully marginalizing over all other 
parameters. As indicated, in the second line BAO data are added to CMB+SN data, in the third 
line also HST data are included. The pivoting redshift is however z p = 0. 



lo c = f2 c /i 2 , 8 = 100 l s /ld, t, n s , log^4, Asz, plus f u when needed (respectively: reduced 
baryon density parameter, reduced CDM density parameter, 100 times the ratio between 
sound horizon at recombination and its angular diameter distance, optical depth due to 
reionization, primeval spectral index, logarithmic fluctuation amplitude with pivoting scale 
0.05 Mpc -1 , SZ template normalization, neutrino fraction as defined below; h is the Hubble 
parameter in units of 100 (km/s)/Mpc). We however kept il^ = 0. 

Our data set includes CMB data from WMAP7 4 from Union2.1 survey [7] (option 
with no systematic errors), WiggleZ [8] and SDSS [9] BAO's data, HST data [10] and 
CMB lensing as provided by CosmoMC. We use different combinations of these data, as 
suitably detailed below. 

In Figure 1 we show la and 2<r contours for the marginalized likelihood on the wq- 
w a plane, when z p = 0, for the sets of data indicated in the frame. In comparison with 
the analogous curves shown in WMAP7 report [12] for f v = 0, our ellipsoids are slightly 
displaced towards more negative w$ and greater w a . The ranges found are closer to the 
Union 2.1 report [7]. 

To gauge the widening of K7o and w a intervals when M u ^ is allowed, we kept the 
same abscissa and ordinate ranges in both sides. The widening is confirmed when fully 
marginalizing in respect to all other parameters, as is shown in Table 1. Let us however 
notice that, when M„ = is required, the inclusion of BAO and/or HST data in top of 
CMB, causes a displacement towards smaller values of the mean wo and an increase of w a . 
These shift - just below la - do not occur (or are much smaller) when M v ^ . 

In Figure 2 we show the likelihood distributions on the f u ~wo and f v -w a planes, 
outlining a progressive delving of w a into the negative domain when f v shifts from to 
0.04 (i.e., when M v shifts from to ~ 0.60 eV). The known result that a greater M v 
is allowed, when w delves in the phantom area, therefore affects w a rather than wq, so 
indicating that, to soften M v limits, it is preferable that w shifts below -1 just when z > . 



4 Provided by the website lambda.gsfc.nasa.gov 
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Figure 3: 2-cr marginalized likelihood contours on the wq — w a plane for different values of z p , 
as indicated in the frame; wq and w a exhibit uncorrelated errors for z p = 0.35 (0.25) when /„ = 
0). Notice the color and linetype inversion between z p = 0.25 and 0.35, in the two Figures. All 
plots refers to CMB+SN+BAO+HST constraint combination. 



3. Results for z p ^ 



Let us then consider the fits when pivoting redshifts z v ^ are considered. In Figure 
3 we overlap the 2-cr contour ellipsoids on the WQ-w a plane, for z p = 0, 0.25, 0.35 and 
0.5, both for f v = (l.h.s.) and ^ (r.h.s.). For the sake of clearity, we consider only 
the full set of observational constraints (CMB+SN+BAO+HST). The ellipsoids exhibit 
a progressive straightening of the symmetry axes and WQ-w a errors become uncorrelated 
when z p = 0.35 or z p — 0.25, in the cases f v — or +~ 0. Such difference between the two 
cases is a result of this analysis. Notice again that the overall ellypso'idal areas are almost 
double on the r.h.s. This is the effect of adding just one extra parameter and confirms 
that a significant correlation exists between the allowed WQ-w a domains and M u , so that 
neglecting the i/-mass option can be badly misleading, when we aim to constrain the DE 
state equation. 
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Figure 4: 2— a marginalized likelihood contours on the wo~w a plane, obtained with different z p and 
translated to a unique redshift. All plots refers to CMB+SN+BAO+HST constraint combination. 
Plots at the l.h.s.(r.h.s.) are translated to z = (to the redshift where wo and w a constraints are 
uncorrelated) . The 2 plots on the first (second) raw are for f v = (allowing f v ^ 0). 

The ellipsoids in Figure 3 apparently distort and migrate through the plot. This is 
because each coefficient pair wo-w a (i.e. WQ Ap -w ai a p ) corresponds to a different straight 
line, when a p varies. 

It is however possible to translate the constraints found at any a p into constraints on 
the plane spanned by the coefficients WQ-w a when z p = 0, or when z p yields uncorrelated 
errors. In the Figures 4 we show the shapes of the ellipsoids after this translation. The two 
Figures on the first (second) line refer to f v = (f u ^ 0). The Figures at the l.h.s. (r.h.s.) 
are the ellipsoids at z = (z = 0.35 or 0.25, where errors on wo and w a are uncorrelated). 

These plots are one of the results of this analysis. They confirm the high reliability of 
the MC algorithm, yielding close results when different parameter combinations are fitted. 
Slight discrepancies however exist: 

(i) In the proximity of the lowest w a values, the limits obtained with different z p exhibit 
some discrepancy. In the f u j^0 case, when z p = 0.50, the major ellipsoidal axis exceed 
other z p by ~ 10 %. 

(ii) The major discrepancy, however, concerns the ellipsoids obtained for f u = and z p = 
(first raw). This is visible in both l.h.s. and r.h.s. plots. In the latter, the wq (w a ) 



- 5 - 



Massless Neutrinos 



Massive Neutrinos 



Wo =t la ± 2a w a ± la ± 2cr ± la ± 2u 



±la±2a f u ±la± 2a 



+o. 

-0. 
+0. 
-0. 
+0. 
0. 

.071°; 
.051°; 



1.13 
1.11 

-1.08 



4.03 



-o 



08+0.34 
10-0.30 
06+0.23 
07-0.22 
05+0.16 
05-0.16 
04+0.14 
04-0.14 
04+0.13 
03-0.14 
05+0.13 
03-0.18 



0.31 
0.32 
0.31 



Ml. 
-0. 
+0. 
-0. 
+0. 



0.321° 
0.281°, 
0.301°, 



40+1.00 

25- 1.32 
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24-1.32 
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-1.10 
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10+0.22 
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56+1.37 
34-1.89 
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37- 1.82 
57+0.40 
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-0.009- 
-0.036- 
-0.009- 
-0.035- 
-0.010- 
-0.036- 
0.009- 
-0.035- 

o n35 +a01 °- 

U.UOO_ 035 _ 

-0.010- 
-0.035- 



0.036; 
0.035; 
0.036 
0.035 



0.035 



0.041 
0.036 
0.040 
0.035 
0.041 
0.036 
0.041 
0.035 
0.040 
0.035 
0.040 
0.035 



Table 2: Mean and fully marginalized limits at la and 2a for different pivoting redshifts. In all 
the cases the data set is CMB+SN+BAO+HST. 




Figure 5: Mean, la and 2a limits on wq at different z p , for massless or massive neutrinos as 
indicated on the frame. In all cases, the data sets is (CMB+SN+BAO+HST). 



range shifts from -1.27/-0.81 (-1.65/1.72) to -1.33/-0.74 (-1.81/1.80) with a unexpected 
area increase by almost ~ 37%. 

The latter finding requires explanation. One could guess that the MC algorithm has 
been optimized for the first and simplest case it had to treat. This guess, admittedly, is 
not supported by any indication pointing at specific options in the operations done. One 
could however deduce from this finding that other wo~w a constraints are probably loose 
and the intervals found are to be reduced by some 10 %. 

4. Discussion 

Linear DE state equations w(a) are expressed by using u>o,a p arid w a ^ ap coefficients (usually 
wq and Wq) which do depend on the pivot a p value selected. Data yield constraints on 
f^o.op and w aAp which turn out to be uncorrelated for a single z p value. We confirm that 
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Figure 6: Envelop of DE state equations (eq. 2). The plot is built by using u> and w a values with 
uncorrelated errors. Any straight line completely inside the black contours is an allowed DE state 
equation. Using other pivoting z p causes almost irrelevant changes (apart of the case z p = with 
fu = 0) . For the sake of comparison, the green dashed lines limit the wq interval allowed when we 
fit data by using varying a p . The data sets is (CMB+SN+BAO+HST). 

the pivoting redshift yielding no correlation is 0.35, if /„ = 0, but we find that it lowers to 
0.25 if a degree of freedom allowing for the i/-mass is open. 

The wo~w a ranges found in the two cases exhibit significant overlaps. If we just 
consider the allowed wq range, when marginalizing in respect to any other parameter, 
we find different trends in the two cases: the u>o range tends to decrease or increase when 
the option f v 7^ is allowed or disregarded. As a consequence, the wq intervals show no 
overlap (at 1 a) when the pivoting redshift exceeds ~ 0.4. This point, as well as the overall 
trends, are shown in Figure 5 and Table 3. 

The set of linear w(a) laws should however be independent of the pivoting redshift. In 
Figure 6 we show the line envelops for the cases f u = (l.h.s.) and 7^ (r.h.s.). In the 
same plots we report also the wq constraints shown in the previous Figure, and the best 
fit wq at any redshift. A hint of the trends found in the evolution of the wo range (and wq 
best-fit value) can be seen also in the average behaviors of the fitting linear laws. 

The fact that the linear laws are the same, indipendently of the z p chosen to perform the 
fit, is confirmed in Figures 4. The plots confirm the good performance of the MC algorithm, 
yielding overlapping results for different linear combinations of the fitting parameters. 
There is however an exception: the fit for the f u = case, performed with z p = 0, yields 
constraints apparently more restrictive (by ~ 37%) than any other z p . 

5. WMAP9 release 

After the completion of this work, nine-year WMAP data have been released, together 
with a number of fits on cosmological parameters [13] (WMAP9). As expected, these data 
yield (slightly) more stringent parameter estimates. 

Let us mention, in particular, the results on wq, w a parameters shown in Fig. 10 of the 
release. The ellipsoids colored in red are to be compared with the blue ellipsoids in Figure 
1 (l.h.s.) of this work. (This is the only significant comparison allowed). 

We had mentioned that our curves are displaced towards significantly more negative 
wq and greater w a values, in respect to WMAP7 [12]. WMAP9 ellipsoids fully confirm 
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such displacement. According to our Figure 1, the 2-a curve in the WQ-w a plane spans 
the intervals -1.5/-0.72 and -1.6/1.6, respectively. The same intervals, in Figure 10 of the 
WAMP9 release are -1.5/-0.84 and -1.2/1.4, respectively. The interval width is then reduced 
from 0.78-3.2 to 0.66-2.6 (for wq and w a , respectively) with an improvement ranging 
between 15 and 19%. The new ellipsoid is completely comprised inside the one shown in 
Figure 1 and almost equally centered with an area reduction by ~ 35 %. In particular the 
shift, in respect to WMAP7 results, due to fresh SNIa data, is fully confirmed. 
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